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Abstract. By treating the electron-ion interaction as perturbation in the first- 
principles Hamiltonian, we have calculated the density response functions of a fluid 
alkali metal to find an interesting charge instability due to anomalous electronic density 
fluctuations occurring at some finite wave vector Q in a dilute fluid phase above the 
liquid-gas critical point. Since |Q| is smaller than the diameter of the Fermi surface, 
this instability necessarily impedes the electric conduction, implying its close relevance 
to the metal-insulator transition in fluid alkali metals. 
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1. Introduction 

The metal-insulator transition in fluid alkali metals such as Rb and Cs has long been 
attracting attention [ll [2], partly because this might be a faithful manifestation of the 
"Mott transition" in its original sense [3] or the one driven by an insufficiently-screened 
Coulomb potential and partly because this accompanies the liquid-gas phase transition, 
the same situation as first considered in mercury by Landau and Zeldovitch [1]. 

In the density-temperature phase diagram of the alkali metals, the liquid-phase 
region is bounded above by the critical point and below by the triple point, as in the 
systems of rare-gas atoms [5]. Above the critical point, there exists the fiuid-phase 
region, an interesting phase in which the metal-insulator transition occurs. In many 
respects, the properties of the dense supercritical fiuid are not very different from those 
of the liquid and a well-known charge instability develops near the transition from liquid 
to solid, as the electronic (ionic) density, n, increases or as the Wigner-Seitz radius of 
electrons (ions) in atomic units, r^, decreases. We can expect, however, to observe a 
different type of charge instability closely related to the metal-insulator transition as 
increases in the fluid phase. 

Conventionally, liquid alkali metals are studied by treating the bare electron-ion 
(pseudo-)potential as perturbation to provide an effective ion-ion interaction via the sea 
of valence electrons described by the three-dimensional electron gas (3DEG) [6]. The 
exchange-correlation effects of the valence electrons play an important role in the dilute 
fluid with large r^. It is fortunate that in developing the flrst-principles calculation 
based on the density functional theory, useful information on 3DEG is now available in 
a wide range of r^. In this paper, we investigate the density response of the fluid alkali 
metal to flnd the charge instability for low densities by use of the standard perturbation 
approach, together with this information on 3DEG. Throughout this paper, we employ 
atomic units. 

2. Response Functions in the First-Principles Hamiltonian 

The alkali metals composed of N electrons and ions are well described by the following 
Hamiltonian H as 

(1) 

Uee = ^ J]Kc(r'j -r^v), 

J2K,ir, - R,,) , (2) 

where {vj} and {Rj} denote, respectiviely, electronic and ionic coordinates, {pj} and 
{Pj} corresponding momenta, m and M the masses, Vee{r), Vai^R) and Vei{r) are the 
bare electron-electron, ion- ion and electron-ion interactions. At long distances Vii{R) 



H = T, + T, + U,, + [/ii + [/ei , 



with 



N 2 
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N p2 

j 

2M' 



Towards the metal-insulator transition in fluid alkali metals 



3 



as well as Ki(^) is represented by a purely Coulombic form, but it deviates from it at 
short-range distances due to the van der Waals attraction and the Born-Mayer repulsion 
between ionic cores for V^i(-R) or due to orthogonality between the valence- and core- 
electron wave functions for Vei{r); the former contribution is weak for alkali metals and 
the latter can be well captured by adopting a suitable (local) pseudopotential \4i(r) = 
Vps{r). In general, the Fourier transforms of these interactions Vai3{q) can be written in 
terms of an effective valence, Za{q), characterizing the long-range Coulomb interaction 
and a short-range part, Uapiq), as 

Vo^M = zmMqHq) + , (3) 

where a and P denote species of the particles, e (electron) or i (ion). Here Zf,{q) = — 1, 
Zi{q) = \yps{q)\/v{q), Ua/siq) = unless a = P = i, and Uu{q) represents the short- 
range interaction between ionic cores. Reflecting the purely Coulombic nature of both 
Vii{R) and Ki(^) at long distances, Zi{q) 1 and Uii{q) converges to a finite value in 
the limit of g = \q\ — > 0. 

Since the first term in the right-hand side of ([3]) is proportional to v{q) = Air/q"^ and 
gets singular at g — 0, it is useful to write the partial number- density response function 
Xa^iq, 1^) in terms of its proper part Ua/si^q, uj), which is irreducible with respect to v{q) 
but not to Uaf3{q), as 

Xapiq, uj) = -Uo^fsiq, Lu) - v{q) ^ Ua-yiq, uj) Z^{q) Zs{q)x5i3{q, uj) , (4) 

7,5=e,i 

which is rewritten as 

/ N rr ^ ^ I ^ ^ "S:^ Ilcy{q,uj)Z^{q)Zs{q)Usf3{q,uj) , . 

X.,{q, u:) = -U^,{q, u) + v{q) 1 + ,(g)n,,(g, ^) ' 

7,(S=e,i 

where the polarization function is defined as Iizz{,q,^) = J2af3 Za{q)Zi3{q)Ilai3{q,uj). 
In the two-component Coulomb system with the effective valences Za{q), the charge 
density is given by the sum of ionic and electronic valence densities, while the (total) 
number density by the sum of ionic and electronic number densities. Then, by use of 
([5D, the charge-density response function xzz{q,(^) = Yj^p ^MZp{q)Xap{q,^), the 
number-density response function XNN{,q,^) = TlapXapiq^^)^ and the cross response 
function XNz{q,^^) = Yliap^P^^)Xai3{.q^^) are, respectively, obtained as 

Xzz{q,^)=-^, — rYFi — ? ^' W 

1 +v{q)Uzz{q,^) 

, ^ n \ I v{q)^Nz{q,^)^zN{q,uj) 
XNN[q,^) = -^NN{q,^) + — ^ — — 1 — ^ — ' (^) 

l+v{q)Ilzz{q,^^) 

llNz{q,uj) 

XNziq, UJ) = --— ^-p^ r , (8) 

l + v{q)Ilzz{q,(^> 

where IlNN{q,^^) = Y^ap na/3(g,u;), YiNziq,^^) = Yap ^/3(9)na/3(q,c^), and IizN{q,i^) 

= Yap Za{q)Ilap{q,Uj). 

Some comments are in order for the so-called g-limits of these response functions. 
In 3DEG or an unperturbed electron system for ([T]), the charge neutrality condition and 
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the compressibility sum rule lead, respectively, to the relations limg^o Xee{q, 0) = and 
limg^o ^ee{q, 0) = n^K with the superscript (0) implying "unperturbed" with respect to 
Uei, where k is the compressibility of 3DEG. These relations reflect a special feature of 
a single-component Coulomb system with a rigid compensated-charge background such 
as 3DEG in which we cannot distinguish the charge-density response function from the 
number-density one. In the two-component Coulomb system, on the other hand, there is 
a difference between them and the difference leads to important modifications on these 
relations. More specifically, the charge neutrality condition leads to 

limxzz(q,0) = limxjvz(q,0) = 0, (9) 

while the compressibility sum rule relates the g-limit of xnn{q,^^) to the isothermal 
compressibility Kt of the total electron- ion system. Then, from ([5]) and (El), it is not 
hard to see that XapiQy^^) satisfies the following relations for arbitrary a and f3 [7]: 

limx«/3(g,0) = ^\imxNN{q,0) = -u^Kt- (10) 

q^O 4 q^O 

Note that the g-limit of XeeiQ-,^) does not vanish but takes a finite value —u'^Kt, 
although limq^ox^e {q,0) = 0. By ffTOl) . we can identify the hquid-gas phase transition 
from a singularity in XaisiQ,^), particularly in the ionic number- density response 
function XaiQ: 0) at g = 0, because Kt diverges at the critical point. Note also that by 
the singularity implying a charge instability can occur only at a finite q. 



3. Approximate Response Functions Based on Effective Ion-Ion Interaction 

Following Ashcroft and Stroud [6j, we can derive an effective ionic Hamiltonian = 
Ti + Uii + Nuq by applying a perturbative method with respect to Uei in ([I]), where uq is 
an energy shift independent of {Pj} and {Rj} (but dependent on n). The interaction 
term Ua can be described by a pairwise sum of effective ion-ion interactions as 

Un = lj2^v.in,-R,,). (11) 

Here the Fourier transform of the effective ion- ion interaction V^i(R) is given by 

VM = Z,{qfv{q)/e{q) + , (12) 

with e{q) the static dielectric function of 3DEG defined as 

e{q) = l + v{q)Iif^{q,Q). (13) 

Equation (fT2|) indicates that the long-range part of the bare ion-ion interaction V^i(q), 
the first term in the right-hand side of ([3]), has been screened by the surrounding valence 
electrons, while the short-range part Ui\{q) remains intact. 

We shall calculate the ionic structure factor S\i{q) in this electron-ion system by 
using the effective ionic Hamiltonian. Since ions can be considered as classical particles. 



-n..(9.°)+ . . L , (15) 
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the classical version of the fluctuation-dissipation theorem relates Sii{q) to the static 
response function Xii(9)0) as 

Xn{q,0) = -nSM/T. (14) 

In the same spirit as in (jlj), we shall introduce Haiq, 0), the proper part irreducible with 
respect to v{q), to write Xii(9)0) as 

0) = -[n,(q, 0)-' + Z,iq)Mq)/eiq) V 

v{q)Z,{qfUy,{q,Or 

l + v{q)[U^\q,0) + Z,{q)m,,iq,0) 

By comparing (fT5|) with ([5]), we find that nee(q,0) = nee''(q,0) and nei(q,0) = 
llie{q,0) = at least in the present approximation. Substituting these equations and 
f fTi|) into ([5]), we obtain 

X.l(9.0) = Xlc(<7,0) = -Z,(q) (l - -ly^ nS„(q)/T. (17) 

From (IT^ . (fT6|) and (ITTI) . the static charge- density, number-density and cross response 
functions are, respectively, given by 

Xzziq, 0) = f 1 - ^) - ^^nSM/T , (18) 



xM,, '^-~^)['-7iqj)-['^ '^^'^ - ) ^^"^^^/^ ' ^^^^ 

0) ^ ^ (l - ^) - ||(l + Z.(.) - ) nSM/T . (20) 

We can easily check that these response functions satisfy the required conditions of ([9]) 
and ([inD with lim^^o 0) = -nSii{0)/T = -n^Kr- 

In what follows, for simplicity, we take Vii(-R) = |-R|~^ and Vci{r) = vji^\r), where 
Vpi^\r) is the Ashcroft's empty-core pseudopotential, given by 



K^Hr) = \ " < ? (21) 





r 


[ -\r\-' 


r 



c 



with Tc being the radius of the ionic core; this simplification leads to Zi{q) = cos(grc) 
and Uii{q) = 47r[l — cos^(grc)]/g^. For e{q), we use an accurate parameterization of the 
diffusion Monte Carlo data of 3DEG at T = by Moroni et al. P]. 
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4. Comparison of Mean-Field Approximation with Monte Carlo Simulation 

We shall now solve the effective ionic Hamiltonian to obtain Sii{q) in the mean-field 
approximation, one of the standard methods in the theory of simple liquids, together 
with Monte Carlo simulations. 

Suppose that the effective ion-ion interaction, Vu{R) = {27i^R)^^ dqVii{q)q sm{qR) 
with R = |R|, can be expressed as the sum of a "reference" part, fo(-R), and a "pertur- 
bation", w{R), by 

V,,{R)=vo{R)+w{R). (22) 

In the mean-field approximation, the free-energy functional F[ni\ of the ionic number 
density ni{R) for the system of interest, characterized by the full potential Vai^R), is 
simply related to that of the reference system -Fo[ni] by 

Flm] = Fo[ni] + J ni{R)w{R - R')ni{R')dRdR . (23) 

Then, we can expand the thermodynamic potential fl[ni] = F[ni] — /ii / ni{R)dR (with 
fii being the chemical potential of ions) up to second order with respect to the ionic 
number-density fluctuations ni(q) as 

fi[ni] = n{n) + ^Y1 SM-'n,{q)ni{-q) + ■■■ . (24) 

Here >S'ii(g) is the ionic structure factor in the mean-fleld approximation, given by 

S,,{q) = [So{qr'+nw{q)/T]-\ (25) 

where So{q) is the structure factor of the reference system and w{q) is the Fourier 
transform of w{R). 

The result of the present mean-fleld approximation depends on how to separate 
the interaction into a reference part vo{R) and a perturbation w{R). A number of 
separations have been proposed so far for the Lennard- Jones potential; among them, 
here we adopt the manner of Weeks, Chandler and Andersen, usually called the WCA 
separation [9]. In this separation, as shown in flgure [H the interaction is split at 
R = Ro, the position of the minimum of Vii(-R), into its purely repulsive and almost 
attractive parts; Vo{R) = Vii(-R) — Eq and w{R) = Eq for R < Rq, while Vo{R) = and 
w{R) = Vii{R) for R > Ro, where Eq is the minimum value of Vn{R) at R = Rq. The 
reference system with the purely repulsive interaction Vq{R) can be effectively treated 
as a hard-sphere fluid where the diameter of the hard sphere is deflned as 



oo 



d= / (l-e-"°(^)/^)di?. (26) 





Using the analytic solution of the Percus-Yevick equation for this hard-sphere fluid 
[Tni[IIl[l2], we can obtain Soiq) and therefore Sniq) by (1251) . 

Figure [2] shows the ionic structure factor Sii{q) calculated in this mean-fleld 
approximation based on the WCA separation (the broken curves) and by Monte Carlo 
simulations (the solid curves) for temperatures and densities corresponding to the 
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conditions of the recent experiments for a liquid Rb [TH] . Note that for all the values of T 
and rg in figure [2] the Monte Carlo results are in good accord with the experimental data, 
which has not been shown here. Although the quantitative disagreement between mean- 
field and Monte-Carlo results can be observed in the vicinity of the triple point (e.g. T 
= 373K, rs = 5.39) and in the vicinity of the critical point (e.g. T = 2173K, = 7.56), 
we can find that the mean-field approximation reproduces at least qualitative features 
of the Monte-Calro results for Sii{q) such as the peak position. In the next section, with 
this knowledge of the mean-field approximation, we will investigate the charge- density 
response of the system in a fluid phase above the critical point. 

5. Charge Instability in a Fluid Phase above the Critical Point 

As discussed in section [21 a charge instability is signaled by a divergence in Xii(9)0) 
= —nS\i{q) /T aX, q = Q = \Q\ with Q being a certain finite wave vector, while the 
compressible instability by that at g = 0. In fact the mean-field approximation based 
on the WCA separation leads us to the result that S\\{q) diverges for low densities as 
shown in figure [3l from which the divergence is seen to occur at g = Q ~ 0.2 for = 14 
(the broken curve) and = 16 (the chain curve), but at g = for = 12.3 (the solid 
curve). Figure m presents the "phase diagram" in the r<j-T plane in which the charge 
instability occurs along the solid curve, while the compressible instability does along 
the broken curve, the former preceding the latter with decreasing temperature in the 
shaded region of > 12.3. 

To obtain a deeper insight into this charge instability in the dilute fluid phase of 
> 12.3, let us consider the charge- density modulation induced by a negative point 
charge put into the origin of the system. This charge- density modulation is given by 
5p{q) = —Xzz{Q,0)v{q) in the linear response theory, which is a total of Spe{q) = 
-Ze{q) Y.a Xeaiq, 0)Z^{q)v{q) and Spi{q) = -Z;{q) J2a Xia{q, 0)Z„(g)w(q) with 6pe{q) 
and Spi{q) being electronic and ionic induced charges, respectively. From ([Hj), ( fT6l) and 
ffTTj) . we can write down 

Spiiq) = f-. nSii{q)/T. (28) 

Figure [5] presents 6pe{q) (the broken curve), Spi{q) (the chain curve) and Sp{q) (the solid 
curve) at two typical points in the r^-T plane marked by (a) the circle with dot and (b) 
the triangle with dot in figure HI we are interested in the charge-density modulations for 
low densities in the former, but for comparison we also show those for high densities in 
the latter which indicate a well-known charge instability related to the transition from 
liquid to solid. In both cases of (a) and (b), the charge neutrality condition 6p{0) = 1 
is satisfied by the fact that the external point charge is completely screened by the 
electrons (|5pe(0)| > |(5pi(0)|), although the ions are highly compressible (|(5pi(0)| ~ 5) 
in (a), while almost incompressible (|5pi(0)| ~ 0) in (b). 
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A couple of notable differences exist between (a) and (b) on this charge instability: 
(1) The "charge-ordering" vector Q as defined by the peak position in Sp{q) is less than 
2pp, the diameter of the Fermi surface of electrons, in (a) but larger than 2pp in (b); (2) 
the charge instability is driven by the electrons, while the ions only screen insufficiently 
the electronic "charge ordering" in (a) because of \dpe{Q)\ > \Spi{Q)\, but it is done by 
the ions with partially-screening electrons in (b) because of \6pi{Q)\ > \5pe{Q)\. 

These differences imply that the charge instability for low densities has a significant 
influence on the electronic transport, while that for high densities does not. It is 
clear that the transition from liquid to solid keeps the system metallic, while it makes 
the translation symmetry broken; this can be understood by the fact that the charge 
instability for high densities is due to the localization of ions with the charge-ordering 
vector Q, which corresponds to the reciprocal lattice vector of the solid with a definite 
direction, and that the Fermi surface of electrons remains almost unchanged since 
Q > 2pp. On the other hand, the charge instability for low densities is characterized 
hy Q < 2pp and then electrons on the Fermi surface connected by Q have strong 
scattering amplitudes with a tendency to localization, leading at least to the decrease in 
the electronic conductivity. Particularly in a fluid or gaseous phase with the translation 
symmetry (in a different manner from incommensurate-charge-density-wave formation 
in a solid), the entire Fermi surface may be seriously affected by Q with randomly 
distributed directions to preserve its rotational invariance, suggesting a new route to 
understanding of the metal-insulator transition. 

6. Conclusion 

We have investigated the static response functions of the fluid alkali metal by applying 
a perturbation method with respect to the electron-ion interaction to the first-principles 
Hamiltonian and found a charge instability in the dilute fluid phase above the liquid- 
gas critical point in the mean-field approximation based on the WCA separation. This 
charge instability is due to electronic density fiuctuations at a finite wave vector Q < 
2pF, implying its intimate connection with the metal-insulator transition in fluid alkali 
metals. 
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Figure 1. Effective ion-ion interaction Vii{R) for — 8.0 and r^ — 2.4 (the solid 
curve), which is divided into the reference part, VQ{q), (the chain curve) and the 
perturbation, w{q), (the broken curve) in the WCA separation. 
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Figure 2. Ionic structure factor 5'ii(<7) in a liquid phase of Rb (rc = 2.4). The soHd 
and broken hues represent Monte-Carlo and mean-field results, respectively. 




Figure 3. Inverse of the ionic structure factor 5ii(<7) in the mean-field approximation 
based on the WCA separation for low densities of = 12.3, 14.0 and 16.0 with the 
core radius of Tc = 2.4. 
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Figure 4. Charge and compressible instabilities in the mean-field approximation based 
on the WCA separation for Tc = 2.4. The structure factor Sii{q) diverges at q = Q 
(with Q beeing a certain finite momentum) along the solid curve, while it does at 
q — along the broken curve in the rg-T plane. 
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Figure 5. Fourier transforms of the electronic, ionic and total charge densities induced 
by a negative point charge at the origin, which are indicated by the broken, chain and 
sohd curves, respectively, (a) for T = 4500[K] and = 15.0 and (b) for T = 2500[K] 
and = 4.0. 



